SUBMODELS OF THE SPECIAL COMPRESSIBLE FLUID ON TWO-DIMENSIONAL SUBALGEBRAS 

1. Introduction 
The differential equation of gas dynamic (EGD) are: 



Cn ■ D^dt+u-\/ 
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Du+-Vp = 0, 

' Dp + pdivu — 0, (1) 
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DS = 0, 



where u is speed, p is density, p is pressure, 5* is entrophy. These variables with 
the special state equation admit 11 parametrical algebra operators Lee Lu. In the 
Cartesian system of coordinates the Lu bases is the following [1, see also 2]: 
C I Xi — dx, X2 — dy , X3 — dz , = tdx + du , ^5 = idy + , 

= tdz + dw, X-! = ydz - zdy + vd^ - wd^, 
Xs = zdx - xdz + wdu - udw, Xq = xdy - ydx + udy - vdu, Xw = dt, 
^ ■ Xii ^ tdt + xdx + ydy + zdz, 

I We considered the special state equation of the kind: 

.O; p^±p-' + FiS), (2) 

where +p'' for 7 > and — for 7 < 0, 7 — , F{S) is the function of entropy. It 
coordinates with the fixed state equation for the fluid at high pressures and high 
temperatures. 

EGD with the equation (2) admit additional operators: 

- stretching X12 = tdt — udu — vdy — wd^ — (7 — '^)pdp — ip9p, 

- carry X13 = dp, 

' where' 7 = 27/(7 — 1), 7 7^ 1- 

Together with Lu L they make up algebra Lee L13. 
' In cyhndrical coordinates (C) x — {x,r,6),u — {U,V,W),y — rcos9,z — 

I rsm9,u = U,v = T^cos^ — Ws'm9,w = ysin6' + Wcos9 the basis of algebra 

Lis is the following: Xi = dx, X2 = cos9dr — sin 6'r~^(9e + Wdy — Vdw), 
X3 = sin6'ar + cos9r-^{d0 + Wdy - Vdw), X^ = tdx + du, 
^ ■ = cos9{tdr - dv) - sm9r-H{de + Wdy ~ (V - rt-^)dw), 

Xe = siii9{tdr + dy) + cos9r'H{de + Wdy - (U - rt-^)dw), X-j = de, 

= '&\^9(rdx-xdr + VdlJ-Vdy)^'CO%9{Wdu-Vdw~xr-^(df) + Wdr-Vdw)), 
^ ; = -co'&9(rdx~xdr^Vdu-\Jdy)+'sm9{Wdv-Vdw~xr-^{dQ^Wdy-Vdw)) 

Ph. 



XiQ = dt, Xix = tdt + xdx + rdr, 

X12 = tdt - Udu - Vdy - Wdw - (7 - '2)pdp - ipdp, Xi3 = dp. 
For algebra L13 all subalgebrus are fisted [4], If 



' 7 = —1,1/3 there are more subalgebrus than for any 7. We shall considers 

two-dimentional subalgebrus from optimum system for L13, appering only when 
7 = —1, 1/3. We shall write out not similar two-dimentional subalgebrus for this 
purprose: 

2.1'. Xi + X2, aXi + Xi3, a(7 - 1) = 0; 
2.2'.Xi2,aXi + Xi3,a^Q; 
2.3'.Xi + Xi2,aXi + 6X5 + X13; 
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2.4'. - Xii + + aXs + Xi3, a 7^ 0; 

2.1".aXi + X12, Xio + Xi3, a ^ 0; 
2.5'.oX7 - 6X11 + X12, Xi + Xi3, a ^ 0; 

2.6'.aX7 + 6X11 + X12, CX4 + Xi3, c2 + (6 + 1)2 ^ 0Va2 + cV 0, 0(7 - 1) = 0; 

2.7'.Xi + aXr + X12, 6X4 + X13, a 7^ 0; (3) 

2.8'.oX7 + X12, 6X4 + Xi3, a 0, 6(7 - 1) = 0; 

2.9'.aX7 - (7 + l)Xn + X12, 6X4 + Xio + X13, 6(7 - 1) = 0, 7 7^ 1; 

2.2".6Xi + aX^ + X12, Xio + X13, h ^ 0; 

2.10'.aX7 + Xio - Xii + X12, 6X1 + Xi3, h ^ 0; 
are for subalgebrus 2.1" and 2.2", 7 -1 7 = 1/3, 
for other subalgebrals 7 = 1 =4>7 = —1, 

here parameters a and 6 which define set of series of dissimilar subalgebras. 
2. Proposal for the coordinates of the equation (2) with the fixed state equation 
The equation (2) coordinates with the fixed state equation [5] 

p = $(p-i)+T/(p-i), (4) 

for same implication functions F{S), f{p^^). Where is describes the 

potential component of pressure, Tf{p~^) is the describes the thermal components 
of pressure, - is the specific volume. The equation (4) describes behavior of 
the real environments, which by their properties approximate firm or fluid frames. 
It is possible with high pressure (about lO^kg/cm^) and high temperature (about 
W'^K). Let us find the meanings F{S), ^(p-^), f{p^^). 

Comparing p (2), (4) and excluding T with the help of thermodynamics (p, S - 
are independent parameters) we receive the identity: 

±pT + F{S) = $(p-i) + {G's - F'sp-^)f{p-% eqno{h) 

where G{S) is determined by the additional experiment. 
After twice differenting by S: 

= -F's- F'^sVfiV) + G'^sfiV), (6) 

where V = p~^. 

1°. Let us Fss 0, then: 

lL^.VfiV) + §^fiV). (7) 
^ss ^ss 

Once again we differentiate on S and receive: (-^-)' = {^^Y f{V)- If 

^ss ss ^ ss 

then, dividing by variables wc have /=const=/o and after the integration wc subtitle 

it in (7). Wc take the contradiction , that p, S are independing parameters. 

Means = 0, i.e. 

^ss 

G'^s = koF^s,Ps = kiF"SS, (8) 

and from (7) follows, that h = -Vf{V) + kaF{V). 

The integration of (8) at fci ^ and the substitution to (5) gives: 



F{S) = kik2e^ + ka, 
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$(p-^) = ±p- + fc3-Tr^- (9) 

Kq — p 



pko - 1 

G{S) = kokikic'^ + UaS + k^, 

where kj are the constants of integration. 

2°. Let us Fss = (it is equivalent A:i = 0). Then F{S) = kiS + ko and from 
(6) we receive (if Ggg ^ 0) -^t- = f{V)=const=fo- Hence from (5) are follows: 

$(p-l) = fcip- Vo + fco - fc2/0 ± 

fip-') = fo, (10) 
G{S) = ^S^ + G,{S) + Go, 

where Go, G'i=const. 

3°. Let us Fss = 0, Gss = 0. Then from (7) follows, that Fs = 0, F{S) = Fq. 
And from (5) we shall receive: 

$(p-i)=±p^ + Fo-Gi/(p-i), 

G{S) = G^{S) + Go, (11) 

where Fq, Gq, G\ are constants. 

Thus, the state equation (2) coordinates with the equation (4), if the functions 
F{S), f{p~^), $(p~^) are represents in one of kinds: (9), (10), (11). 

3. Calculation of invariants 

For construction the submodels of the special compressible fluid necessary to 
calculate invariants of the subalgebrus. [2, see also 1]. 

The algorithm of the calculation invariants consist in the following: 

1. Wc select the system of coordinates, in witch calculates the invariants. If 
the subalgebra contains the operator Xr of the rotation, it is convenient to choose 
cylindrical coordinates. If the operator of the rotation id not resent are convenient 
the cartesian coordinates. 

2. We write out the operators of subalgebra in convenient system of coordinates 
from the list (3). 

3. We enter the function h, which depends from 9 variables {t,x,u, p,p) as 
required invariants. 

4. The function h is invariants of the subalgebra L =< Yi,Y2 > Only when any 
operators Y of subalgebra, working on invariants function, to annul it. Namely, 
Y ■ h = 0,Y € L. We shall work by the operator Yi of the basis subalgebra L on 
invariant function. In result we received the linear homogeneous equation with the 
partial derivatives of the 1-st order. For this equation we write the characteristic 
equation, the system of the ordinary differential equation [6]. Let us assume, that 
there is an obviously complete set functionally independent invariant (integrals) 
l''{t,x,u, p,p), k = 1..8. 
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4. We write down the second operator of the basis through received invariants 
by the rule: 

Y2=ed.^=C~dj. (12) 

5. We shall work by the stayed operator Y2 on invariant function h{I^). We 
received the linear homogeneous equation with the partial derivatives of the 1-st 
order. Wc write down for it the equation of the characteristic. We find a complete 
set of invariants. 

6. We pass to initial variable. 

The received invariants (3) are shown into the table (sec the appendix). 
Example: 

As an example we shall consider the subalgebra 2.7' from (3): 
Yi = Xi + aX^ + X12 = ade + tdt - Udu - Vdy - Wdw + pdp - pdp, 
Y2 = bX4 + Xi3 = btd^ + bdu + dp. 

Let us the invariant function h(t, x, u, p,p), x — (x, r, 6), u = {U, V, W), 

satisfying to the equations Yx-h = 0, Y2-h = 0. 

The second equation looks like hthx + hhu + /ip = 0. 

Let us write down the equation of the characteristics: 

dx _dU _dp _dr _de _dV _dW _dp _ dt 

We find integrals, which form a complete set functionally independent invariants : 
t;p;W;V;e;r;Ui = U -xt'^wi = p-x{bt)-^. 

We shall receive hix = 0, when written down the second equation through in- 
variants by a rule (12). Means h = hi{t, r, 9. V, W, p,pi,Ui). 

The first equation looks like in new invariants for the known equations variable 

a/ii0+t/iit+(-[/4<ct~^)/iicri-y/iiy-W/iivr+p/iip+(-p-l-a;(6t)~^)/iipi=O; 

Having written down the characteristic equation and having calculated integrals, 
we receive a complete set functionally independent invariants, which in initial vari- 
able look like : 

r;e- a\n\t\; Ut - x; Vt; Wt; pt~^;pt - xb~^. (13) 

4. Invariant submodels of the second rank 

Two-dimcntional subalgebra has 5 invariants. The invariant decision is exist if 
all required functions are defined from expressions for invariants. These invariants 
are nominated by new functions from others invariants for this purpose. Others 
invariants necessarily will be the functions of independent variables [2]. All un- 
known function are defined from the received equality. Thus, the representation of 
the invariant decision turns out, which is substituted to EGD. The system of the 
equations turn out which connect only invariants and new invariant functions as a 
result of the substitution under the theorem of the representation invariant variety 
[6]. The equations for invariants refers an invariant's submodel. 

We shall write down an invariant submodel for the considered example. 

We shall make equality from invariants (13): 9 — a\n.\t\ = 9i, 
Ut-x = Ui{r,ei),Vt = Viir,9i),Wt = Wiir,9i), pf-^ = pi{r,9i),pt - xb'^ = 
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The representation of the invariant decision is defined from these equality: U = 
{Ui+x)t-^; : V = Vit-^; : W = Wit''^ ; p = pit; : p = pit'^ + x{bt)-'^ 

The representation of the invariants decision for S can be received from the state 
equation: p = ±p~^ + S = t~^{x{b~^) + Si), where Si =pi±p~^ replaced the 
state equation in the invariant submodel. 

The statement in EGD results to the invariant submodel: 

Di = {Wir-^ - a)dg, + Vidr, 

DiUi = -{pih)-\ 
DiVi+pirPi'^ =Wi\-^ + Vi, 
DiWi+pie,{pir)-^ = Wi- ViWir-\ (14) 
Dipi+pi{Vir + r-^Wie,) = -pi{2 + Vir-^) 
DiSi = -Uib'K 

Any invariant submodel can be resulted to the one of the two initial types by 
the choice invariants [7]: 

- Evolutionary (time-f is the invariant of the subalgebra) 

D = dt+ U2ds, 

Du2 + hpi~^pis = ai, 
Dv2 = a2, 

Dw2 = as, (15) 
Dpi + piU2s = ai 
DSi = as, 
6>0; 

- Stationary 

D = U2da;^ +V2dyi, 
Du2 + bipi~'^Pi^^ = ai, 

Dv2+b2Pl~'^Plyi = 02, 

Dw2 = as, (16) 

Dpi + Pl{U2xi + V2yi) = a4, 

DSi = as, 
6i > 0, : 62 > 0; 

here Oj, b, bi is the coefficients of the initial types. 

From examining subalgebras (3) is received two submodels of the evolutionary type, 
and from others subalgebras is received ten submodels of the stationary type. 

The canonical types of invariant submodels are tabulated (see appendix), where: 

- 1-st column is the number of subalgebra, 
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- 2-nd column is the basic system of the coordinates in which considered the 
EGD, 

- 3-rd column is the initial type: S is the stationary type , E is the evolutionary 
type, 

in 4-th column are given invariants, 

in 5-th column the factors of the initial type are written down. 

The (14) implements to stationary initial type by replacement 
r = xi. : 6 — alii\t\, : U2 = Vi, : ?,'2 = ixi)^^Wi + a, : W2 = Ui with factors: ai = 
U2 + xi{v2 + a), : (12 = {v2-a)(l-2u2{xi)^'^), : 03 = l-(p6)^^,a4 = -pi{u2Xi^ + 
2), : as = (1 + W2)b-^ + Si, : 61 = 1, : 62 = 

The example of reduction subalgcbra 2.9' to the canonical type. 

The operators of the subalgebra are those: 

Yi = aXr - 2Xn + X12, 

Y2 = 6.T4 + Xio + Xi3, &(7) = 0, : 77^ - 1. 
Invariants from independent variables look like: Xi = {x — b2^^t^)r~^ , 
yi = 6 + a2~^ln|r|. The representation of the invariant decision enters the name 
through the new invariant function: V = Vir^,:W ~ Wir^, : p = pir~^,p = 
Pir2 -\-t, : U = Uir2 + bt, where Vi,Wi, p,pi,Ui depends on X\,y\. 

The representation for entropy is defined from the state equation 5 = S\r^ + 1, 
where Si = pi ± pi~^. 

The substitution to the EGD results the following invariant submodel: 

Di = {Ui - xiVi)d,, + (Wi + a2-^Vi)dy„ 

DiUi + pr^ = -b- 2-^ViUi, 

DiVi + pi-\piy,a2-^ -Pi^.xi) = W^-pi{2pi)-^ - 2-^Vh 
DiWi+pi-^Piy,=WiVi, 

Dipi + piiUi,, - Vi^.xi + ViyM'^r)-^ + Wiy,) = -3Vipi2-\ 
DiSi = -1 - ViSi2-'^. 
The new invariant speeds are entered on expression for Di: 

Ui - XiVi = U2, : a2~'^Vi + Wi = V2, Wi - 2a-^Vi - xiUi2a-'^, with which we 
receive replacement: X2 = xi^—ayi, : j/2 = yi-|-2~^aln U3 = 2xiU2—av2, ■ Vs = 
{2xi)~^au2 + V2- From which follows the system (16), where: 

ai = -2.T2& - Vi{x2Ui + 2x2{Ui - xiVi) - aWi + ViX2{p)-^) + {2x2 - a){W? - 
pi(2pi)-i - 2-1T/2) ^ 2{Ui - xiVi)^ 

a2 = (a(2x2)-i + l)(W^i^ -pi(2pi)-i - 2-^1^) - 06(2x2)"^ - aViUi{Ax2)-'^ + 
{2x2)-^a{Ui - 0:1^1)^1 - 2-iax2-2 + WiVi, 

03 = WiVi-2a-\Wi^-pi{2p)-^-2-'^Vi^)+2x2a-^{-b-2-'^ViUi)+2a-\Ui- 
xiVi)Ui, 

tti = -^piVi, 

as = -1-2-11/1^1, 

bi = {2x2^2-^a^f + 1 + 4x22, 62 = 3^ + 1, pi = ±pi-i + Si. 

5. Invariant submodel of the third rank 

The expressions for invariants arc define the speed and the pressure, but it is 
impossible to define density (see table, appendix) for subalgebra 2.1" from optimum 
system (3) at o = 0. It is possible to build a regular partially invariant submodel 
in this case. 

Let is give the definition to the regular partially invariant decisions generally. 
Let us for algebra H are present 7i, .., 7^- invariants from independent variables 
and are present Ji,.., J;- invariants from dependent variables. If from invariants 



7 



Ji,..,J/ are defined all dependent variables, it is possible to build an invariant 
submodel of the rank k, nominating invariants Jj by functions from (Ji, .., 1^), i.e. 

Jj = Jjih,..,h),j = l,..,l. (17) 

If it is impossible all dependent variables of invariants Jj, to define, then (17) 
gives the representation of the regular partially invariant decision of the rank k, 
and defect a,, which is equal to number not determined independent variables, i.e. 
a = m — I, where m is the number of dependent variables. 

For subalgebra 2.1" the rank equal 3, the defect is equal 1. 

Let us consider the subalgebra 2.1" more in detail. 

The operators of basis are those: 

Yi=dt + dp, 

^2 = tdt — udu — vdy — wdw + 'ipdp + pdp. 

Invariants from independent variables are : x,y,z. From others invariant, spec- 
ified in the table, we receive the representation of the regular partially invariant 
decision.. 

u = p^ui{x,y,z), : p = t + p^pi{x,y, z), : p = p{t,x,y,z). (18) 
The substitution to the EGD, gives: 

1 1 2 1 1 

-3'"i(Pt + p'^uvVp) + pi[{uvV)ui + Vpi] + -p~^pvVp = 0, (19) 

2 1 2 

Pt + ■^p'^uv'Vp + p^divui = 0. (20) 

shall received the representation of the decision for entrophy from the state 
equation S = t + pa Si, where 5*1 = pi — 1. 

The substitution to the equation DS = 0, gives: 

1 2 1 

-Sip-s {pt + p-^uvVp) + uvVSi + 1 = 0. (21) 
From (20) and (21) are follows: 

^ = -9^-^±^^+3clivui. (22) 
p Si 

Then it is possible to find pt from the (21): 

ipt = 3p-i[-divui + 2Si'^{ui-VSi + l)]=p-iB{x). (23) 

Replacing pi on 5i + 1 and substituting (23), (22) to the (19) we receive: 

-•Vp = [-^ui{l + ui-VSi) - VSi - {ui-V)ui]-\-=A{x). (24) 
p Si Si + 1 
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By the substitution (24) to (22), we exclude p: 

(^J^ - 3)(«rV5i + 1) + ^^^{u,.V){S, + \u,^) = 0. 

Equating the mixed derivative functions Inp from (23), (24), we receive VB = 
^BA, rotA = 0. Prom the last equality follows, that A = V(p and V(31ni? — i^) = 

O^SlnB - if = 0=>B = esv. Prom (24) follows p = b{t)e'^. Then from (21) we 
receive 6' = b^. 

The integration gives b= (|)^, where constant of the integration is made by zero 
with the help of carries on t and on p, admitted by EGD. 

So, density as p = t^pi {x, y, z), i.e. is the representation of the invariant decision 
for one-dimensional subalgebra Y2. 

Thus, the is a reduction of the partially invariant decision to the invariants: 

(Ul • V)W1 + -pi = Ml, 

Ui ■ V/9i + pidivui = — 3pi, (25) 
Ml • 5i = -Si, 

where Si = pi — pi^ , S = tSi. 
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Appendix 

2.1'. C.S.: D, Type: S; 

Invariants: y, z, tv, tw, tu — x + In \t \,pt~^, tS — xa~^ + \n\t\; 
Submodel (16): ai = U2,a2 = V2,a3 = 1 — {api)~^,a4 = — 2pi,a5 = 51-1- 
a"i(l -'W2),bi = 62 = 1; 
2.2'. C.S.: D, Type: S; 

Invariants: y, z, tv, tw, tu — x, pt~^, tS — xa^^; 

Submodel (16): ai = U2,a2 = 1^2, as = —{api)~^,a4 = — 2pi,a5 = 5i -|- 
a~^W2,bi = 62 = 1; 

2.3'. C.S.: D, Type: S; 
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Invariants: x — ayb ^ — In | t \,z,tu — ab ^tv — l,tw,ab ^{tu — ayb ^ — 1) + 
tv -y,pt-'^,tS -ya-^; 

Submodel (16): oi = U2 + a{b'^p{)~^ + l,a2 = ^2^1,03 = ab~^{u2 + 1) — 

04 = — 3pi, 05 = Si — ba~'^W2 + a~^U2, bi = a?b~'^ + 1, 62 = 1; 
2.4'. C.S.: D, Typo: E; 

Invariants: t, z^^{y — atx), z~^{v — ax — atu — sw), z~^[{a'^t'^ + .s^)(u — ax) + 
atu + sw], z~^{(a^t'^ + s'^)w + s{v — ax) — atsu)pz, z^^{S — x); 

Submodel (15): ai = -2(1 + a^t^ + S^)-'^[u2{a'^t'^ - stU2) + W2{tu2 - s)] = 
ai (pi ) - 1 - spi (pi ) " \ a2 = (1 + a^t^ + ) - ^ [2ai;2 - 1*2 aH"^ - 2asw2 + (w2 +U2)i2aH + 
SU2) - W2V2 + SV2U2 + W2U2 + s(m2)^] + ai(pi)~^ - spi(pi)~\a3 = (pi)~^(pi(l + 
aH'^ + ats) + {I + a^t"^ + s'^)~'^[2aH'^W2 - (102- su2)H-2sv2 + as'^W2] + (^2)^,04 = 
-2p((w2 - SM2)(1 + aH'^ + s^)-\ as = [w2(S'i - s) +V2- U2{aH + Sis)]{l + aH^ + 
s^)-\b={l + aH^ + s^),Pi=pr^ + St; 

2.5'. C.S.: D, Type: E; 

Invariants: t,e + a\a\t |, r-^(ay + VF), r-^C/, ar"^(Wa - V),pr, r"^(S' - x); 

Submodel (15): ai = -a(pi)-i(5i + [pi)-^) + {a{a'^ + 1))-Hw'2' - a'u2' + 
2u2W2, a2 = -{pi)~^ - V2{a{a'^ + l))~^(aM2 - W2), ^4 = -pi(a(a^ + l))~^(a2W2 - 
W2), (I5 = -V2, 6 = + 1; 



